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1-006 ae tt 9 OTS il YRTHRTL. SRE SATHRANDOM. MAR: 1 . (1) 
-TITLE MTHSRANDOM - coered | Purpose Random Number Generator 
-IDENT /1-006/ 3 File: MTHRANDOM.MAR 


MPARBARAAAAAALAAAAAAAAL ALLA A RESALE ASA RARE ES ERE REESE SRSA REE E ERECT EERE 


:* COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT EORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


ED UNDER A LICENSE AND MAY BE USED AND COPIED 
THE RMS THE 
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:* TRANSFERRED. * 
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® 
® + 
® ~ 
® ® 
® ® 
® * 


3* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


CORPORATION. 


DIGITAL ASSUMES NO RESPONS 
SOFTWARE ON EQUIPMENT WHIC 


ITY FOR THE USE OR RELIABILITY OF ITS 


BIL 
IS NOT SUPPLIED BY DIGITAL. 


I 
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++ 
FACILITY: Mathematics Library 


ABSTRACT: 


This is a general random number generator. It is 
of the mul spe teat ive congruential type, and hence 
is fast, although prone to certain classes of 
non-random sequences. 


ENVIRONMENT: any access mode, normally user mode 
modular, AST reentrant procedure 


AUTHOR: Peter F. Conklin, CREATION DATE: 28-Dec-77 
MODIFIED BY: 


Peter F. Conklin, 1-Feb-78: VERSION 01 
Original creation. 
Change to LIBS rather than MTHS. 

Change back to MTH$. Some more comments. JMT 4-Jan-78 
4 Correct documentation on usage. : 
1-005 = Put version number in standard format (three digits in 

edit number field). JBS 16-NOV-78 
1-006 - Add ‘*_"' to the PSECT directive. JBS 22-DEC-78 
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00000004 
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«SBTTL DECLARATIONS 


; INCLUDE FILES: 

; NONE 

; MACROS: 

: NONE 

; EQUATED SYMBOLS: 


SEED=4 :; Only argument is the seed 
OWN STORAGE: 


NONE 


+ PSECT DECLARATIONS: 
.PSECT _MTHSCODE,PIC,SHR,NOWRT,LONG, EXE 
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4 ag -SBTTL MTHSRANDOM = CALLable function 
: FUNCTIONAL DESCRIPTION: 
RANDOM - General Purpose Pseudo Random Number Generator 


Call the function MTHSRANDOM (seed) to obtain the 

next pseudo-random number. The seed S$ updated by 

the function automatically as a side effec™. The 

result is a floating po nk number that is uniformly 
distributed in the range 0.0 inclusive to 1.0 exclusive. 
There are no costc tes tons on the seed, although 

it should be initialized to different values on 
separate runs. 


The algorithm used is to update the seed as: 

seed = 69069. * seed + 1 (mod 2**32) 
and then to convert the seed to floating point. 
Note, because the result is never 1.0, a simple 
way to get a uniform random integer selector is 
to multiply by the number of cases. For example 
if a uniform choice among 5 situations is to be 
nade, then the following FORTRAN statement will 
work: 

GO TO (1,2,3,4,5) 1+#IFIX(5.*MTHSRANDOM( seed) ) 
Note that the explicit IFIX is necessary before 


adding 1 in order to avoid a possible a lg 
during the normalization after the floating add. 


SOoOoOoCSoooCoooosooo 
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This is a general random number generator. It is 
of the multiplicative congruential type, and hence 
is fast, although prone to certain classes of 
non-random sequences. This non-random behavior 
typically arises when considering triples of 
numbers generated by this method. 


For more information on congruential generators, 
see: 


Random Number Generation (pp. 1192-1197) 
by G. Marsaglia 


in: Encyclopedia of Computer Science 
edite by masnene Ralston 
Petrocelli (New York, 1976) 
CALLING SEQUENCE: 


result.wf.v = MTHSRANDOM (seed.mlu.r) 
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INPUT PARAMETERS: 


SEED. mlu.r Longword seed is modified each call 


IMPLICIT INPUTS: 


NONE 


OUTPUT PARAMETERS: 


NONE 


IMPLICIT OUTPUTS: 


NONE 


FUNCTION VALUE 


Returns in RO a single-precision floating point value between 
0.0 inclusive and 1.0 exclusive. 


SIDE EFFECTS: 


NONE 


Page 
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ALLable function “-SEP-19 


4 MTHRTL. SRC JMTHRANDOM.MAR; 1 


ENTRY MTHSRANDOM,0O ino registers save, clear IV 


If this were to be placed as an inline expansion, then 
GAUL SEED ,#69069,41,R0 should replace the next two 
geyruct ions because this would prevent the possibility 
of integer overflow trapping. 


MULL2 #69069,aSEED(AP) 


supdate seed with multiplier 
INCL @SEED(AP) 


sincrement seed to protect 
3; against strange seeds 


The next instructions convert the seed from unsigned integer 
to floating point in the range 0.0 to 1.0 exclusive. 


EXTZV #8,#24,a@SEED(AP) ,RO :Get the most significant bits 
; of the se d in the range 


CVTLF RO,RO Convert to floating without 
; rounding. The result is 
3 ow beat and in the range 
3 0.0 .. (2.0%#24)-1.0 


: If this were to be placed as an inline expansion, then 
3; MULF #*x00003480,R0 could replace the next two instructions. 
BEQL 10$ sIf zero, already correct 
SUBW = #2487, RO SDIVE #*F2,0%224 
; the result is :.ow in the 
; range 0.0 .. 1.0 exclusive 
10$: RET 
- END 


e working set tare A 900 pages 
1893 bytes (4 pages) of virtual memory were used to buffer the intermediate code. 
there were 10 pages of symbol table space mh ieceses to hold 2 non-local and 1 Local symbols. 
source Lines were read in Pass 1, producin object records in Pass 2. 


pages of virtual memory were used to define 8s macros. 


Macro Library name 


~-5255$DUA28: CSYSLIBISTARLET.MLB; 2 0 


0 GETS were required to define 0 macros. 
There were no errors, warnings or information messages. 
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Paeoemsceaconeerecmneed 

! Psect synopsis ! 

foe new cemmam cmon eee t$ 
PSECT name Allocation PSECT No. Attributes 
- ABS, 0000000 ( ot ( -) NOPIC USR CON ABS’ LCL 
_MTHSCODE taetdt ( 3$:} 8 ( 9.) PIC USR CON REL LCL 

dem woaaweoecenwoee ramen een ena $ 
H Performance indicators ; 

Phase Page faults CPU Time Elapsed Time 
Initialization 2 : 00:0 0.08 :00:00.55 
yomnene processing 1 . :00: “2 :00: +°08 
Syabol table sort : :00: 3 :00: 0:69 
Symbol table output :00. 0: if 219 
Psect synopsis suteut :00: “Oe 0:00: 8°08 
Cross-reference output :00: 09 Se Be 
Assembler run totals 291 0:00:01.5 0:00:08. 
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Page 
MTHRTL.S . (&) 


NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SHR EXE RD NOWRT NOVEC LONG 


MACRO/ENABLE =SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHRANDOM/OBJ=0BJ$:MTHRANDOM MSRC$:MTHRANDOM/UPDATE=(ENHS$ : MTHRANDOM) 


ORPORATION 
OPRIETARY 
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